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ABSTRACT 

Galaxy clusters at high redshift present a superb opportunity to study galaxy evolution, with large galaxy sam- 
ples at fixed distances. Several lines of evidence point towards large formation redshifts and passive evolution for 
cluster ellipticals. At redshifts larger than z ~ 0.8, this picture rests on data from only a handful of rich clusters. 
We have observed four potential clusters to sample sparser environments in a critical redshift range, 0.8 < z < 1 .2. 
We compare the photometric evolution and radial density profiles in our clusters with more massive clusters stud- 
ied in a similar manner. We also highlight the overdense ERO (I — K > 5) population, significantly redder than 
the old elliptical population at the cluster redshift. We discuss the implications of our clusters at high redshift, 
and possible differences in environment from richer clusters at similar redshifts selected through X-ray or optical 
techniques. 

Subject headings: clusters: observations — galaxies: evolution — galaxies: formation — infrared: galaxies — 
radio: galaxies 



> 

o 
o 

(N 

\o ■ 
o : 
o . 
o ■ 

Or 

o ■ 
a: 



1. INTRODUCTION 

High redshift galaxy clusters present luminous signposts for 
the formation and evolution of massive galaxies. While the 
study of field galaxies at very large redshifts (z > 3) has recently 
progressed (e.g. Steidel et al. 1999), clusters continue to be our 
best probe of dense environments and large mass scales at early 
epochs. Clusters with luminous AGNs at their centers offer a 
different perspective on galaxy interactions and activity than 
clusters selected either on the basis of optical over-densities 
or spatially extended X-ray emission. The AGN-selected clus- 
ters tend to be poorer than the optically and X-ray selected rich 
clusters that dominate most catalogs, and thus probe a different 
range of environments (Lilly & Prestage 1987). Such a sparse 
environment may also be related to the formation and fueling 
of the luminous AGN (Heckman et al. 1986). 

The existence of massive collapsed objects, particularly rich 
galaxy clusters, at high redshifts provides a challenge for the- 
ories of cosmic structure formation (e.g. Press & Schechter 
1974). The z~l range represents a critical point at which seri- 
ous constraints can be placed on the formation and evolution- 
ary history of early-type galaxies in clusters. Beyond z~l in 
cosmologically flat CDM models, the amount of merging oc- 
curring within the prior ~1 Gyr was thought to be sufficiently 
large so as to dramatically inflate the locus of early-type galaxy 
colors. Recent modelling has shown however that the locus of 
such red galaxies may remain quite tight until beyond z ~ 1.5 
(Kauffmann & Chariot 1998). However these models still have 
difficulty predicting the correct slope of the red sequence. 

The emerging evolutionary picture from the study of a few 
high redshift galaxy clusters (e.g. Ellis et al. 1997, Stan- 
ford et al. 1998, Kodama et al. 1998) is one of smooth and 
steady spectro-photometric evolution for early-type galaxies 
over more than half the age of the universe - implying large 
formation redshifts and primarily passive evolution. The aver- 
age E+S0 galaxy color at a particular luminosity evolves grad- 
ually, while both the slope and scatter of galaxy colors around 
the mean change very little. The small scatter implies that, even 
approaching z^l, we have not reached the time at which a sig- 
nificant fraction of the stars in early-type galaxies were young. 



The fact that the color-magnitude slope does not change with 
redshift strongly favors a scenario in which this slope arises 
from a mass-metallicity correlation (Stanford et al. 1998). Dif- 
ferent burst ages and epochs would result in a changing slope 
with redshift. 

The evolution of luminous AGN may be tied to the evolution 
of their parent environments, particularly in clusters. The peak 
in the QSO luminosity evolution is thought to have occurred 
between l<z<3 (e.g. Boyle & Terlevich 1998). This was the 
age of quasars. In the present epoch, powerful AGN are far 
less numerous and appear not to inhabit very overdense regions. 
However, prior to z^0.5 conditions were different, many radio- 
loud AGN existed in overdense or cluster environments (e.g. 
Hall et al. 1998, Yates et al. 1989, Yee and Green 1984; Hill 
& Lilly 1991; Ellingson, Yee, & Green 1991). If mergers and 
interactions are indeed the key to fueling luminous AGN, then 
the temperatures and velocity dispersions in later-epoch clusters 
(z ~ 0) may be too large for efficient mass transport. The strong 
QSO evolution may be naturally explained as a result of cluster 
virialization and evolution - at higher redshift, the lower galaxy 
velocities lead to higher interaction probabilities (e.g. Fabian & 
Crawford 1990). A luminous AGN present at the cluster center 
may be an indicator of rapid dynamical evolution of the cluster 
and the constituent galaxies. 

Although a few recent studies have addressed cluster ellip- 
tical evolution at redshifts ~1 (Stanford et al. 1998, Stanford 
et al. 1997, Rosati et al. 1999, Tanaka et al. 2000, Kajisawa et 
al. 2000), there is a marked paucity of data available for clus- 
ters in this critical redshift regime (z~l), especially those with 
less rich environments where mergers are expected to be more 
efficient. We have identified a sample of potential high redshift 
clusters from a complete sample of radio galaxies (Chapman 
et al., in preparation), the Molonglo Reference Catalog 1 Jy 
sample (McCarthy et al. 1996; Kapahi et al. 1998). We have 
selected all sources identified with galaxies with measured red- 
shifts in the range 0.8 < z < 1.5, using visible and /Ts-band 
images obtained with the 2.5 m du Pont telescope at Las Cam- 
panas. Since old elliptical cluster galaxies have very red colors, 
/T-band galaxy selection removes a large source of blue inter- 
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loping field galaxies, highlighting potential cluster members. 
In addition, cluster members do not suffer as greatly from k- 
corrections in the Zf-band as in the optical (Stanford et al. 1998) 
and TT-magnitude/color diagrams therefore present a compara- 
ble study of cluster evolution across a range of redshifts. 

In this paper we present deep multiband images of four ob- 
viously overdense radio galaxy (RG) fields from our sample at 
z > 0.8. Confirming the cluster status will ultimately require 
spectroscopic measurements of redshift for many cluster mem- 
bers. With the present data, we rely on Color-Magnitude dia- 
gram (CMD) analyses to reveal the characteristic red sequence 
of elliptical cluster members, in consort with the measured ra- 
dio galaxy redshift, and in one case (MRC 0959-263) a few 
cluster member redshifts lying fortuitously along the slit (sec- 
tions 3.1,3.2). A radial density analysis is also used to assess the 
cluster richness over the field population (section 3.3). We also 
highlight the overdense ERO population, significantly redder 
than the old elliptical population at the cluster redshift (section 
3.4). In section 4, we discuss the implications of our clusters 
at high redshift, And possible differences in environment from 
richer clusters at similar redshifts selected though X-ray or op- 
tical techniques. 

2. THE SAMPLE, OBSERVATIONS AND REDUCTION 

The four radio galaxy fields in this study were chosen from a 
large sample of K- and V-band images, because of an apparent 
overdensity of red galaxies in the K image when compared with 
V . We also required that the (V-K) vs K CMD revealed a linear 
sequence of red galaxies including the radio galaxy as further 
evidence for the cluster environment. This was not meant to be 
systematic for the sample, but merely to choose a few potential 
clusters at z ~ 1 for deeper followup observations. 

2.1. Imaging 

The optical and near-IR images on which our present analysis 
is based were obtained using the du Pont 100 inch telescope at 
Las Campanas, Chile, using CCD cameras in the optical. The 
J and H-band images were obtained with the CIRSI camera 
(Beckett et al. 1998 - lkx lk pixel Hawaii array), while the Ks- 
band images were taken with IRCAM (Persson et al. 1992 - 
NICMOS3 256 2 array). The details of the optical observations 
and their reduction can be found in McCarthy et al. (1996). The 
near-IR observations were taken with a 30" dither pattern, and 
reduced in a standard approach for IR imaging. 

To facilitate the analysis of the datasets, the various images of 
the four cluster fields were resampled and rotated to common 
reference frames using the IRAF tasks GEOMAP and GEO- 
TRAN. The fields are centered on the initial estimate of the 
cluster positions from the red galaxy spatial distribution. The 
final images have the same pixel size as the 0.265 arcsec/pixel 
CCD images. 

The source extraction software SExtractor (Bertin & Arnouts 
1996) was used to make an initial object detection pass at >3er 
with the final catalog at >4a in the TTs-band images, while mea- 
suring flux parameters for each individual passband in the sep- 
arate images. Magnitudes and colors were measured using a 2 
arcsec diameter aperture, roughly twice the FWHM of the im- 
ages. The depths reached in each filter vary somewhat field to 
field, but are typically close to 4a detection limits of Ks=20.5, 
#=21.4, 7=21.4, 7=24.0, V=25.0. The actual limits for each 
field are presented in Table 1. Figs. 1 and 2 present the V and 
TT-band images of the four clusters, while Fig. 3 depicts a false 
color image of the richest cluster in our sample, MRC 1022- 



299. For visualization, an I — K color-coded representation of 
each field is presented in Fig. 4, with red squares highlighting 
probable cluster members (those lying within the presumed lo- 
cus of early-type galaxies as discussed in the next section). 

The Galactic E(B-V) toward the four fields is small (« 
0.04), and we have not applied an extinction correction to our 
photometry. We estimate the quadrature sum of the systematic 
lcr errors in the photometry (zeropoints and PSF matching) to 
be ^0.05 for the IR colors, and ^0.04 for the optical-^ colors. 

2.2. CTIO 4m spectroscopy 

The RC spectrograph on the CTIO 4m telescope was used 
to obtain optical spectra for the four radio galaxies, although 
no clear redshift was measured for MRC 0527-255. For 
MRC 0959-263, two other galaxies falling on the slit had red- 
shifts within 5% of the radio galaxy, and are likely cluster mem- 
bers. Details of the observations and reductions can be found 
in McCarthy et al. (1992), while the spectra will be presented 
in McCarthy et al. in preparation. 

3. ANALYSIS AND RESULTS 
3.1. Color-magnitude diagram analysis 

The color-magnitude diagram (CMD) forms our main diag- 
nostic for identifying red galaxy overdensities, and we present 
a brief overview of the technique before discussing each cluster 
individually. The CMD for a galaxy cluster typically reveals 
a tight red sequence of massive elliptical galaxies, implying a 
mass-metallicity relation for early-type cluster galaxies (Stan- 
ford et al. 1998, Gladders et al. 1998). This argument can be 
inverted to search for the presense of a cluster via the red ellip- 
tical sequence. The fact that studies to date have implied large 
formation redshifts and primarily passive evolution for the clus- 
ter ellipticals allows the red sequence to be used as an estimate 
of the cluster redshift, even at redshifts out to z ~ 1 . We have 
used this approach to characterize the environments of the four 
radio galaxies under investigation here. In the cases where the 
radio galaxy redshift is known, we have, along with zcmd, two 
independent measures of the cluster redshift to compare. 

Figure 5 shows the CMDs for the RG fields in observed- 
frame colors. For the two higher redshift clusters, J- and H- 
band images were obtained since the near-IR colors begin to 
show a significant deviation from the field population. Two ra- 
dial cuts have been made on the object catalog centered on the 
radio galaxy, with solid squares representing objects within 35 
arcsec of the radio galaxy, crosses from 35" to 60" and circles 
for objects greater than 60" . The solid diagonal lines represent 
4cr detection thresholds in color. A color-magnitude sequence 
of red galaxies, which includes the radio galaxy, is visible in 
Fig. 5 for all clusters, with many of the red objects lying within 
the central 35 arcsec. Our position-color representations of the 
cluster fields (Fig. 4) depict the red sequence objects with red 
squares, while blue triangles depict bluer objects lying below 
the red sequence. 

Plotted in all panels of Figure 5 (dashed lines) are esti- 
mates of the no-evolution color-magnitude locus for early- 
type galaxies. The dashed lines are derived from the dataset de- 
scribed in Stanford et al. (1998) which uses multiband photom- 
etry of early-type galaxies covering the central ~1 Mpc of the 
Coma cluster (Eisenhardt et al. 1997, Stanford et al. 1998). The 
dashed lines therefore represent the colors that Coma galaxies 
would appear to have if the cluster could be placed at the red- 
shifts of our radio galaxies and observed through the VRIJHKs 
filters we used. The dashed lines were interpolated from the 
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values presented in Stanford et al. (1998), where they have a 
cluster very near the redshift of each cluster in our sample to 
Z = 1. For MRC0527-255, we made a very small extrapola- 
tion to the no-evolution Coma line from the Lynx field clusters 
at z = 1-26 presented in Rosati et al. (1998), and Stanford et 
al. (1997). 

Our procedure was to fit a line to all objects within 18 < 
K < 20, taking a color cut of /- K ± 0.4 around the Coma no- 
evolution line derived for the radio galaxy redshift. For this 
line we then determined the equivalent Coma no-evolution red- 
shift. In fitting the red sequence in this manner, we noticed 
that measuring the slope and scatter of these objects is plagued 
with systematic errors, since the values changed considerably 
depending on the I-K region used to fit over, with no objec- 
tive criterion to select the I-K fitting region. For the one RG 
without a redshift, MRC 0527-255, the dashed line represents a 
fit of the no-evolution Coma colors to the brightest red galax- 
ies including the RG, within the central 35 arcsec, satisfying 
the color cut of 3.6 < I-K < 4.4. The results are presented 
in Table 1, and compared to the measured redshifts of the ra- 
dio galaxy in the cluster. The uncertainty in the derived red- 
shift results from neglecting to correct for Galactic extinction, 
from the systematic zero-point uncertainties in the Coma data 
(which we take as 0.075 mag - Stanford et al. 1998), the pho- 
tometry, and the choice of objects used for the fit, giving rise to 
a ~ 0.10-0.15 mag error for the individual clusters. The I-K 
red sequence shows a much tighter relation than V - K, likely 
because the extinction is less and the effects of any recent star 
formation are less pronounced at the longer wavelength. 

The case for MRC 0527-255 is less clear than the other three. 
The I-K red sequence containing the radio galaxy indicates a 
redshift of z ~ 1.3, which would explain the difficulty in mea- 
suring a redshift for the radio galaxy with CTIO 4m blue spec- 
trum. Spectra of the brighter cluster members will be required 
to verify whether this indeed represents such a high-z cluster. 

Certainly for the three radio galaxies with measured red- 
shifts, the location of the expected locus of Coma colors is re- 
markably consistent with the colors of the red excess galaxies 
in the field. MRC 0959-263 and MRC 1 139-285 do show small 
color offsets from those expected from the RG redshift (see sec- 
tion 4). The color difference is measured by the difference at 
K = 19 between the Coma no-evolution line at the RG redshift 
and our red sequence fit line. 

3.2. Individual Clusters 

In this section, we present details of the individual radio 
galaxy (RG) fields in order of increasing redshift. 

3.2.1. MRC 0959-263 

The K=\l.l radio galaxy lies near the bright end of the red 
sequence of galaxies, but a K=ll.O elliptical lying 20.7" to the 
northwest may actually be the brightest cluster member (BCM) 
galaxy. The RG has a redshift of z=0.68, but the extrapolated 
redshift from the red sequence is higher at z=0.75. The no- 
evolution I-K color difference is ^0.15±0.10 magnitudes. 

3.2.2. MRC 1139-285 

This radio galaxy appears to be a likely candidate for the 
BCM, situated as the brightest galaxy on the red sequence. An- 
other bright red galaxy lying next to the RG on the CMD dia- 
gram is actually a barely resolved pair lying 15" to the north- 
west. Although the RG has a redshift of z=0.85, the cluster 
sequence shows very similar colors to MRC 0959-263, with an 



estimated redshift of z=0.77. This is not very significant against 
the no-evolution I-K color difference of ^0.1 8 ±0.11 magni- 
tudes. The cluster members are much more spatially dispersed 
than in the other three fields, with less obvious central overden- 
sity surrounding the RG. 

3.2.3. MRC 1022-299 

A false color image of MRC 1022-299 (Fig. 3, plate 1) shows 
the dramatic overdensity of red galaxies, the most obvious over- 
density from our sample of four radio galaxy fields. This is 
also apparent in the clear pileup of galaxies along the red se- 
quence in the color magnitude diagram (Fig. 5). The red se- 
quence, with an estimated redshift z = 0.95, clearly includes the 
RG. The cluster redshift, estimated from the red sequence, is 
slightly higher than the spectroscopic redshift of the RG itself 
(z = 0.93). The radio galaxy is not the BCM, with the actual 
BCM on the red sequence (ZTs=17.8) lying 38" from the RG. 

3.2.4. MRC 0527-255 

Although no redshift was measured for for the radio galaxy 
MRC 0527-255, the -redshift relation for radio galaxies (e.g. 
Lilly 1989, McCarthy 1993) and the 8" apeture K magnitude 
of 17.2 implies z ~ 1-2, a higher redshift than the other radio 
galaxies in the present sample. However, the K-z relation has 
significant scatter, (dz~0.3 RMS at z = 1.2), and it would not 
be surprising to find the RG and its associated over dense envi- 
ronment at a somewhat lower redshift. 

The radio galaxy is amongst the brightest members of a loose 
sequence of galaxies with very red colors, where the estimated 
sequence redshift would be z = 1 .28. Comparison with the two 
clusters at similar redshift discovered by Rosati et al. (1999) 
and Stanford et al. (1997) show very similar colors in all bands, 
lending plausibility to this field representing a high redshift 
(z > 1) cluster. 

In Fig. 6 we plot the measured flux converted to AB- 
magnitudes in all known bands of the six -brightest objects ly- 
ing on the red sequence, within the central 35" radius of the RG. 
The TT-band flux for all galaxies has been normalized to that of 
the RG. Overplotted are Coleman, Wu & Weedman (1980) el- 
liptical galaxy models redshifted to z = 0.8, 1.3 and 2.0. The 
templates are at least suggestive that if these objects are old 
cluster members, they lie at redshifts greater than z = 1, and are 
consistent with a cluster lying at z = 1.3. 

3.3. Radial Density 

Having identified apparent red sequences from the CMD in 
all four candidate clusters, we can explore the density of the 
fields as a function of radius from the radio galaxy. The galaxy 
density was calculated in successive annuli 5" wide for each 
cluster and for the average of all four clusters together (Fig. 
7a). Figure 7a clearly shows a significant excess above the 
field galaxy limit at our survey depth, Ks=20.5 (solid horizon- 
tal line with lcr errors due to variation in counts across the 
field as dashed lines - Cowie et al. 1994,1996 and Saracco et 
al. 1999) for the four clusters averaged together. The average 
integrated counts over the central 25" radius are > 6a above 
the field population, estimated conservatively as the field count 
plus lcr. However the signal in the core region for any cluster 
on its own is still convincing. 

By taking an I-K color cut on the objects of 1 magnitude 
bluer than the red sequence in each case (Fig. 7b), the signifi- 
cance of the overdensity signal becomes larger, as we have ef- 
fectively removed a portion of the contaminating field galaxy 
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population, (peaking at 65 arcsec" 2 w 6a above the field) How- 
ever, such a cut also removes all the blue galaxy population 
residing in the cluster which was bright enough to be detected 
at Zfs-band. Cowie et al. (1996) present the counts of field ob- 
jects of the Hawaii deep surveys with various I—K color cuts. 
In Figure 7b the field levels are estimated for each cluster using 
the Ks=20.5 count with an I—K color cut one magnitude bluer 
than the red sequence value at Ks=\9 - the central value of our 
fitting procedure described in section 3.1. The I—K field lim- 
its vary for the clusters from ~ 6/arcmin 2 (MRC 0959-263) to 
~ 3/arcmin 2 (MRC 0527-255). The average integrated counts 
of the respective I—K field limits over the central 25" radius is 
a 7.5ct detected overdensity. 

There appear to be two spatial components to the overden- 
sities in all clusters (Fig. 7). The first is an overdensity within 
~30" of the RGs relative to our own data at >30" (a near-field 
excess). For all objects with Ks< 20.5 this represents a 4.5cr in- 
tegrated excess averaged over the four clusters. However, even 
excluding this near-field excess, the integrated counts for all 
clusters at radii >30" relative to random-field counts from the 
literature (Cowie et al. 1996, Saracco et al. 1999) represent a 
4.5cr detection. This far-field excess appears to extend to >80" , 
although our data are sparse at these radial distances. The am- 
plitude of the near-field excess is insufficient to explain the total 
excess Zfs-band counts, and the far-field excess is likely real. 
The overdensity in our four fields corresponds to an approxi- 
mate Abell class of 1.5±0.5, where the red sequence extrapo- 
lated redshift is used in all cases to measure the angular scaling. 

3.4. EROs 

Radio galaxy and quasar fields have shown apparent over- 
abundances of EROs (e.g. Liu et al. 2000, Dey, Spinrad & Dick- 
inson 1995, McCarthy et al. 1992). The most natural explana- 
tion is that these might represent evolved early type galaxies in 
a cluster or group environment. Without morphological or sub- 
mm information, however, it is difficult to know if any might 
represent dusty obscured starforming objects (Dey et al. 1999, 
Cimatti et al. 1999). Although all those objects with I-K > 4 
would be considered anomanously red if discovered in the field, 
in high-z cluster fields they lie close to or along the red galaxy 
sequence, and their apparent overdensity in our four clusters 
is large, especially in the two higher-z fields. However, the ob- 
jects with I—K ;> 5 are clearly too red to form part of the typical 
mass/metallicity related sequence of cluster members. 

We define two samples of EROs in our RG fields: those with 
I-K > 5, and those objects which are significantly above the 
red sequence (I-K > 4.3 for the lower-z clusters, I-K > 4.7 
for the higher-z clusters). Although many such objects are de- 
tected, very few of these objects have K < 19 and a direct 
comparison with the random field ERO(/? - K > 6) population 
of Thompson et al. (1999) is difficult. In Table 1, we list 
the ERO(I-K > 5) and ERO(I-K >red seq. ) numbers for 
each cluster to limits of K < 20.0 (S/N> 5a) and K < 20.7 
(S/N> 3.5er). As some of the candidate EROs appear near the 
edges of the field, where the exposure time is less than in the 
central regions, they do not meet our detection criterion of 5a 
and 3.5a for the two respective samples, and we do not include 
these in our catalog. For ERO(I-K > 5), we find 10 and 18 
objects in the K < 20.0 and K < 20.7 samples respectively, 
while for ERO(I-K >red seq. ), we find 14 (K < 20.0) and 
30 (K < 20.7) objects. The number of EROs is similar in the 
four clusters. As the optical magnitudes are very faint in many 
cases (below our nominal Aa detection threshold), all the ERO 



candidates with Ks< 20 were verified to be present by visual in- 
spection of the multi-filter data. Several of the K < 20.7 sample 
were only detected at ZTs-band. 

A field survey using HST-NICMOS data (Yan et al. 2000) 
derives surface densities of 0.31 ±0.14 arcmin" 2 and 0.63±0.2 
arcmin" 2 , with R-H ;> 5 and H < 20.0, H < 20.5 respectively. 
Yan et al. (2000) liken their R-H > 5 to R-K > 6, which has 
been shown to be roughly equivalent to our I-K > 5 ERO cri- 
terion (Liu et al. 2000). Comparing directly our ERO(I-K > 5) 
samples to Yan et al. (2000), we find surface densities of 
0.78 ± 0.25 arcmin" 2 (2.5 x the field) and 1.41 ± 0.34 arcmin" 2 
(2.2 x the field). 

The EROs are not generally within the ~ 35" near-field over- 
density region, and tend towards the outer regions of our ob- 
served fields (with the possible exception of MRC 0527-255). 
Our candidate EROs are represented as black circles in the 
position-color plots (Fig. 4). This may be indicative of the for- 
mation processes that lead to the extremely red colors if the 
EROs are related to the high-z cluster environment. A small 
percentage of the already overdense old and red cluster galax- 
ies might have enough dust to further redden their colors signifi- 
cantly above the mass/metallicity sequence. Alternatively, there 
may be an overdensity of dusty ULIRG-type galaxies formed in 
the dense environment of the radio galaxy. Their location in the 
outer field excess would be consistent with objects undergoing 
some phase of merger (e.g. van Dokkum et al. 1999). Sensitive 
sub-mm observations should eventually be able to distinguish 
between these two possibilities, since the dust mass in the lat- 
ter case would be substantially larger than for slightly dusty red 
cluster ellipticals. 

In addition to the above EROs, we also found 4 objects which 
originally appeared to be EROs from the I-K color, but were 
found to be bright at V-band. Three of these objects appear in 
the highest-z field, MRC 0527-255. These objects appear simi- 
lar to the FROG objects (Faint Infrared-Excess Field Galaxies) 
discovered by Moustakas et al. (1999), and may represent ob- 
jects with a nuclear starburst or AGN. 

4. DISCUSSION 
4. 1 . Comparison with other clusters 

The magnitude, spatial, and color distributions of the excess 
galaxy population in these RG fields are all consistent with the 
excess being produced by overdensities of galaxies at the RG 
redshifts. Roughly speaking, the amplitude of the near-field ex- 
cess corresponds to Abell richness class ~1 clusters, and the 
far-field excess to Abell richness ~1.5. This is consistent with 
RGs often being located in large-scale galaxy overdensities and 
occasionally in smaller-scale groups within them. However, the 
richness measurements should be interpreted cautiously, as one 
galaxy at z=0 may typically correspond to several galaxies at 
z~l due to mergers in the interim (Steinmetz & White 1997, 
Stanford et al. 1997). 

The objects within 60 arcsec of the radio galaxies at K < 20 
do not have colors which are significantly different from the 
no-evolution Coma colors transformed to the redshift of the 
RG. Because the uncertainty in the zeropoint of the Coma no- 
evolution prediction is similar to the amount of expected color 
evolution in J-K and H-K, we would not find strong evi- 
dence for passive evolution of the cluster galaxies even without 
our uncertainties in fitting the red sequence cluster redshift. 

The lack of redshift identifications for cluster members, and 
the large potential for intervening structure at these relatively 
high redshifts, makes it difficult to estimate color evolution. 
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We have shown that at least for MRC 1022-299, the red se- 
quence estimated redshift is very close to the RG redshift itself. 
For the case of the more dispersed cluster, MRC 1 139-285, and 
the nearest cluster, MRC 0959-263, the red sequence estimated 
redshift is somewhat different from that of the RG (fe=0.08). 
Although the radio galaxy redshift is not necessarily the best 
redshift for the cluster (the RG is not clearly the BCM galaxy 
in any cluster), a velocity dispersion corresponding to dz=0.08 
would be relatively large for clusters at our RG redshifts (Del- 
torn et al. 1997, van Dokkum et al. 1999, Rosati et al. 1999). 
It would therefore be surprising if the cluster redshift were sys- 
tematically offset from that of the radio galaxy. In addition, 
the photometry and Coma zeropoint errors added in quadra- 
ture (3cr ~ 0.2 mag) are similar to the measured color offsets. 
Therefore we are reluctant to associate such minor color offsets 
from the Coma no-evolution colors with excess or diminished 
cluster evolution. MRC 0527-255 has only the red sequence es- 
timated redshift to work from, and until spectroscopic redshifts 
are obtained, nothing further can be gathered about the color 
evolution. 

The colors of most of the objects lying near the red sequence 
are broadly consistent with the predictions of Bruzual & Char- 
lot (1997) elliptical galaxy models for z = 0.76,0.78,0.95 and 
1.28 respectively, with the following properties: passive evo- 
lution of a 1 Gyr burst stellar population with solar metallic- 
ity formed at z = 5 for h = 0.65, qo = 0.1. For our most dis- 
tant possible cluster (MRC 0527; z = 1.28), such a model pre- 
dicts V-K = 7.0, I-K = 4.3 and J-K = 1.9 for a galaxy 
age of 3.25 Gyr. It is noteworthy that the same "standard" 
passive-evolution model was also found to provide a reason- 
able fit to the average optical-IR colors of early-type galaxies 
in the large sample of rich optical and X-ray selected clusters 
at 0.3 < z < 0.9 (Stanford et al. 1998). Ages less than ~3 Gyr 
for the stellar populations could be accommodated if the mean 
metallicity were greater than solar. Also, we have made no at- 
tempt to correct the colors for extinction due to dust internal to 
the galaxies. If this is large, then the true colors of the cluster 
galaxies would be bluer than our measured values, resulting in 
a better fit to models with z/ < 3. 

The apparent scatter and slope of the red sequence can be 
used as a measure of the relative mix of stellar populations 
present in these early epoch cluster ellipticals. Our lack of mor- 
phology selection from high spatial resolution (ie: HST) im- 
ages makes color scatter comparisons difficult, since we have 
inadvertently removed all bluer galaxies from consideration as 
cluster members. Stanford et al. (1998) noted no significant in- 
crease in red sequence scatter to z ~ 0.9 using HST-selected 
E/S0 galaxies in clusters. Beyond z= 1, the mix of stellar 
populations is expected to increase the scatter and slope (Ko- 
dama et al. 1998), leaving a clear excess of red galaxies but a 
much looser "sequence". This appears to be the case in two 
recent analyses of AGN-selected clusters: 3c 324 at z = 1.206 
(Kajisawa et al. 2000) and B2 1335+28 at z= 1.11 (Tanaka 
et al. 1999). However, the z ~ 1.3 clusters discovered in the 
Lynx IR field survey, RXJ0848. 9+4452 and CIG J0848+4453 
(Rosati et al. 1999), still show a relatively tight red sequence. 
Our apparent red sequence scatter and slope out to z ~ 0.9 ap- 
pear similar to the eye as that found by Stanford et al. (1997), 
but we do not attempt to quantify the slope or scatter as de- 
scribed in section 3.1. MRC 0527-255, possibly at z ~ 1.3, does 
show a rather loose red sequence similar to 3c 324, although the 
pile up of red objects at fainter magnitudes is fraught with pho- 
tometric errors, and it is not clear which objects to include in 



the sequence. 

4.2. Connections with the radio galaxy 

These four clusters are clearly amongst the highest den- 
sity environments in our total sample of high-z radio galaxies 
(Chapman & McCarthy, in preparation). The radio power and 
source structure varies significantly through the four RGs, even 
after accounting for redshift differences (Table 1). However the 
spectral index (0:843/408) is similar for the four objects, placing 
them all slightly shallow of the steep spectrum source defini- 
tion (a > 0.9) from Kapahi et al. (1998). The relatively steep 
spectrum is typical for distant RGs (McCarthy 1993). The RG 
is not clearly the BCM galaxy in any of the clusters (although 
it is one of the two or three obviously brightest sequence mem- 
bers in MRC 0527 and MRC 1139). This may be expected if 
the formation and evolution of the radio galaxy depends for in- 
stance only on a few trigger merger events early in the cluster 
formation period. The radio galaxy properties would then only 
loosely relate to the local environment. 

Although early suggestions that nuclear properties and radio 
source structure are linked to the clustering environment (Heck- 
man et al. 1986), we find no obvious evidence for environmen- 
tal differences in our cluster fields over those selected through 
X-ray or optical techniques. Higher merger rates might be ex- 
pected to associate with high-z RG cluster environments, espe- 
cially those less rich clusters with smaller velocity dispersions 
(e.g. Barnes 1999). However, direct studies of merger fractions 
are difficult without higher resolution imagery, and blue frac- 
tions of clusters at high-z are difficult to quantify due to the 
strong variation in the blue field population in the foreground 
(Ellis 1997). Our clusters do all show a compact environment 
present in the central ^35" which is somewhat distinct from 
richer X-ray selected clusters at similar redshifts (van Dokkum 
et al. 1999, Rosati et al. 1999). This may indicate a locally 
higher merger environment in the cluster center than typical 
(by contrast, van Dokkum et al. 1998 found higher merger frac- 
tions in the outer radii of the clusters through HST imaging of 
close pairs). It remains to be seen whether there is any true cor- 
relation of cluster and radio properties in a large flux limited 
sample. 

5. CONCLUSIONS 

The observations presented here provide strong evidence 
supporting the cluster identification for the very red galaxy 
overdensities in the environments of these four radio galaxies 
environments. The density of all objects within 25"radius of 
the radio galaxy field centers is ~ 6 x higher than the field aver- 
age. Excluding the bluer objects, the overdensity peaks at ~ 7.5 
times the field count subjected to similar /- K color cuts. 

The quiescent behavior of early type galaxies in clusters to 
z ~ 1 appears to represent the evolutionary history of these RG- 
selected clusters fairly well. The smooth and steady evolution 
picture has been demonstrated with only a handful of relatively 
rich clusters selected through the X-ray and optical morphol- 
ogy thus far, and our four clusters add significant evidence to 
this picture. 

No clear connection exists in our data between the radio 
galaxy properties and the clustering environment, except for the 
local peak in the overdensity surrounding the radio galaxy. 
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FIG. 1 . — Cluster field images comparing the near-IR to optical. Top row: MRC 0527-255 K, MRC 0527-255 V; field size 95 arcsec on a side. 
Bottom row: MRC 0959-263 K, MRC 0959-263 /; field size 1 10 arcsec on a side. The radio galaxy is marked in the ^s-band images. The other cluster 
members with confirmed redshifts in MRC 0959-263 are circled. In all fields in this paper, north is up, east is left. 




FIG. 2.— Cluster field images comparing the near-IR to optical continued. Top row: MRC 1022-299 AT, MRC 1022-299 V; field size 90 arcsec on a 
side. Bottom row: MRC 1 139-285 K, MRC 1 139-285 V; field size 90 arcsec on a side. The radio galaxy is marked in the ATs-band images. 
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FIG. 3. — V,I,K-band false color image of the richest cluster in our sample (120 arcsec on a side), MRC 1022-299. The excess red galaxy population 
clearly stands out in an already overdense field. 
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MRC 0527-255 MRC 0959-263 




FIG. 4. — The 4 cluster position-color diagrams showing the density of red sequence and ERO objects. The angular extents of the plots are similar 
to the images in figures 1&2. Only objects with a > 4cr detection at A^s-band are depicted. Top row: MRC 0527-255 and MRC 0959-263 Bottom 
row: MRC 1022-299 and MRC 1139-285. Solid red squares represent objects lying on the red sequence, solid blue triangles represent galaxies with 
I — K > 1.2 mag. bluer than the red sequence, and open blue Mangles with I—K > 2 mag. bluer than the red sequence. EROs with I — K at least 1 
mag. redder than the respective elliptical sequence are shown with black circles (at least I—K > 4.3). 
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FIG. 5. — The cluster CMDs at 2 radial cuts showing the red sequence emerging. In order of increasing redshift; Top row: MRC 0959-263; 
MRC 1 139-285, Bottom row: MRC 1022-299; MRC 0527-255. Radial cuts have been made on the objects to provide an idea of the color distributions 
with distance from the radio galaxy. Solid squares represent objects within 35 arcsec of the radio galaxy, crosses from 35" to 60" and circles from 
60" to 90" . Dashed lines depict the locus of Coma cluster colors transformed to the redshift of the RG-selected cluster. The solid diagonal lines 
represent 4a detection thresholds in color, with magnitude errors added in quadrature. Note that for MRC 0959-263, the Gunn r-band image covers a 
somewhat smaller area than the / and K images, and many of the galaxies identified in the I,K fields are not present in r. 
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FIG. 6. — MRC 0527-255: A spectral energy plot of the six if -band brightest galaxies lying along the assumed red sequence for the cluster. The lines 
represent an elliptical galaxy redshifted and k-corrected to z = 0.8 (dotted), z = 1 .0 (dashed), and z = 1 .28 (long dashed). Although evolutionary effects 
will make this model SED less red, dust effects would considerably redden the profile, providing a better fit to the data at wavelengths shortward of 
1 micron. These six red objects provide some indication that the cluster likely does lie at z > 1 . 
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FIG. 7. — Fig. 7a (left panel). Projected radial distribution of all galaxies brighter than K s = 20.5 (2" aperture diameter) within the average of the 
4 clusters (filled circles), with error bars are calculated for the number of objects in each bin as per Gehrels (1986). The individual clusters are also 
plotted, slightly offset in radius for clarity, with MRC 0527 (triangles), MRC 0959 (squares), MRC 1022 (pentagons) and MRC 1139 (hexagons). A 
clear signal is seen as a significant excess above the field galaxy limit at our &>-band survey depth (solid horizontal line, la error bars as dashed lines). 
Fig. 7b (right panel). By making a color cut on the objects of 1 magnitude bluer than the I—K red sequence in each case, the strength of the overdensity 
signal becomes even larger, as we have cut out a portion of the contaminating field galaxy population. The I—K color-cut limits for the field in each 
case are taken from Cowie et al. 1996 (horizontal dashed lines), and vary from ~ 6/arcmin (MRC 0959-263) to ~ 3/arcmin 2 (MRC 0527-255). Again, 
the error bars are calculated for the number of bins as per Gehrels (1986). 
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Table 1 

Properties of the Radio-Galaxy Clusters 



Property 


MRC 0527-255 


MRC 0959-263 


MRC 1022-299 


MRC 1139-285 


Ks 4a limit 


20.5 


20.3 


20.5 


20.4 


H 4a limit 


21.2 


n/a 


21.3 


n/a 


J 4a limit 


21.0 


n/a 


21.4 


n/a 


/ 4a limit 


24.6 


24.4 


24.5 


24.5 


V 4a limit 


25.5 


24.6" 


25.5 


25.2 


RG redshift 


n/a 


0.68 ± 0.005 


0.93 ± 0.005 


0.85 ±0.005 


red sequence redshift 


1.28 ±0.03 


0.75 ± 0.02 


0.94 ±0.02 


0.77 ± 0.02 


Abell class 


1.5 


1.5 


2 


1 


No. red galaxies, R<10" 


6 


5 


3 


2 


No. red galaxies, R<35 


22 


13 


16 


9 


No. red galaxies, R<60" 


32 


31 


38 


20 


ERO,_k >5 to #=20.0(20.7) 


3(6) 


3(5) 


2(3) 


2(6) 


ERO> rec | set |+i mug 


4(10) 


3(7) 


4(5) 


3(8) 


S mM Hz(mJy) 


1540 


1720 


1120 


6810 


Q 843/408 


0.88 


0.82 


0.84 


0.85 


radio structural type 


partly resolved 


triple source 


partly resolved 


double source 



REFERENCES. — a) This object has gunn-r magnitude rather than V 



